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Abstract. We investigated the MeV properties of 173 unidentified or only tentatively identified EGRET sources 
listed in the third EGRET catalogue, by analyzing the simultaneously collected COMPTEL MeV data for each 
individual source. The sources can generally be divided into 4 groups. In this paper we focus on one of these, a 
group of 22 EGRET sources for which we can provide additional constraining information: their spectral extrapo- 
lations from the energy range above 100 MeV towards lower energies overshoot the fluxes or upper limits derived 
simultaneously at MeV energies. This means that for these sources a spectral turnover/break between 1 MeV and 
100 MeV is required. At least two of these sources, but most likely the majority of this sample, have the maxima of 
their gamma-ray luminosities in this energy band. The sources have rather soft EGRET spectra (average photon 
index = 2.72lg;?f), and seem to spatially cluster in the inner Galaxy. Variability analyses revealed 11 out of the 22 
sources to be significantly variable. Object classes proposed as possible counterparts for the unidentified EGRET 
sources are discussed in the light of these additional constraints. 
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1. Introduction 

One of the biggest mysteries left by the Compton Gamma- 
Ray Observatory (CGRO, 1991-2000) is that a large num- 
ber of 7-ray sources detected by the different CGRO ex- 
periments, in particular EGRET, still remain unidentfied. 
The EGRET experiment measured 7-rays above 30 MeV, 
most sensitively above 100 MeV. Out of the 271 sources 
listed in the third EGRET catalogue ijHartman et all 
1999), 170 are unidentified and 27 are only tentatively 
identified. Several classes of objects have been proposed 
as possible counterparts for those unidentified EGRET 
sources. Sources located at high galactic latitudes and 
being time variable are believed to be active galactic 
nuclei (AGN), in particular blazars. Sources, located at 
lower galactic latitudes, being steady and having low 
7-ray fluxe s, are found to coi ncide spatially with the 
Gould Belt (|Gehrels et all 120001. Some other low-latitude 
sources show positional correlations with supernova rem- 
nants (SNRs) and OB associations (e.g. iRomero et all 
1999). Steady sources with hard 7-ray spectra seem to 
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be good candidates for young 7-ray pul sars with ages of 
less than 10 6 years l|Zhang et all EoOOaj) . Several sources 
(mainly located |b| < 10° 1 might indicate a new class of 
7-ray emitting objects ((Torres et all l200llh because they 
do not coincide with any potential counterpart objects. 

The COMPTEL experiment aboard CGRO is sen- 
sitive to 7-ray photons between 0.75 and 30 MeV, 
thereby covering the softer 7-ray band adjacent to the 
EGRET one. Apart from transient 7-ray bursts, uniden- 
tified 7-ray sources and AGN are the majority of the 
COMP TEL source detecti o ns. T he first COMPTEL cat- 
alogue l|Schonfelder et all l2000)l lists 10 AGN and 9 
unidentified 7-ray sources; the sum of the rest (radio pul- 
sars, stellar black-hole candidates, SNRs, and 7-ray line 
sources) is about 12. Since COMPTEL and EGRET were 
mounted parallel on CGRO and both had a large field of 
view (the COMPTEL one being larger than the EGRET 
one), COMPTEL and EGRET observed simultaneously 
the same sky region. 

To gain further knowledge on the unidentified EGRET 
sources, and to probe their nature, we analyzed the 
contemporaneous COMPTEL data on the unidentified 
EGRET sources to supplement the EGRET results. In 
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this paper, we report the discovery of a subgroup of the 
unidentified EGRET sources whose 7-ray spectra are con- 
strained by the MeV data: their spectral energy distribu- 
tions have at least an MeV break but most likely an MeV 
peak. The paper is organized as follows: in Section 2 we 
briefly describe the COMPTEL instrument, the applied 
data analysis methods and the observational concept of 
CGRO, in Section 3 we present the analysis results and 
discuss them in Section 4. In Section 5, we finally present 
the conclusions. 

2. Instrument, Data Analysis and CGRO 
Observations 

The Compton telescope COMPTEL (0.75-30 MeV) had 
an energy-dependent energy and angular resolution of 5% 
- 8% (FWHM) and 1.7°- 4.4°(FWHM), respectively. Rs 
field of view is circular and covers ~1 steradian. Imaging 
in its large field of view is possible with a location accuracy 
(flux dependen t) of the order of 1 ° - 3°. For details on the 
experiment see ISchonfelder et al.l i|l993h . 

Skymaps and source parameters, like detection signif- 
icances, fluxes, and flux errors, can be obtained via the 
maximum likelihood method, which is implemented in the 
standard COMPTEL data analysis package. The detec- 
tion significance is derived from the quantity -21nA, where 
A is the ratio of the likelihood L (background) to the like- 
lihood Li (source + background). The quantity -21nA has 
a Xi distribut ion, if only the flux at a given source position 
is estimated ijde Boer et all Il99^) . The detection signifi- 
cance can be conservatively calculated by the ratio of flux 
to flux error. This approach is adopted in this paper for 
estimating the source detection significances. 

In fitting the fluxes of the relevant EGRET sources, 
nearby prominent COMPTEL sources are taken into ac- 
count by fitting simultaneously their fluxes. An estimate 
for the instrumental background of COMPTEL is derived 
by using the standard filter t echniq ue in the COMPTEL 
data space (jBloemen et all ll994jl . The celestial back- 
ground components, galactic and extra-galactic diffuse 7- 
ray radiation, also have been taken into account by model 
fitting. In the presented analyses we applied instrumental 
point spread functions assuming an E~ 2 power-law shape 
for the source spectra. 

CGRO observations were organized in so-called 
'Mission Phases' and 'Viewing Periods (VPs)'. A 'Mission 
Phase' covers typically 1 year of data and contains many 
VPs, which typically last for 1 to 2 weeks each. For each 
EGRET source of interest we analyzed the simultane- 
ously collected COMPTEL data. To allow combining of 
COMPTEL and EGRET results, we analyzed the MeV 
data for periods for which the EGRET spectral index was 
estimated. Many EGRET source results are published in 
the sum of the CGRO Mission Phases 1 to 4, noted as 
'P1234' in the third EGRET catalogue. For this time pe- 
riod we generated COMPTEL all-sky data, and derived 
the flux results by fitting sources at the relevant positions. 
As noted above, the fit was performed by including 1) 



models for the cosmic diffuse radiations, 2) the strongest 
MeV source (the Crab), and 3) neighboring COMPTEL 
sources on a case by case basis. 

To compare our COMPTEL source fluxes with the 
EGRET ones, we plotted the best-fit power-law shapes 
to the EGRET source spectra with la errors and extrap- 
olated these below 100 MeV into the COMPTEL band. 
Systematic errors of 10% are included in these EGRET 
spectra. The COMPTEL flux values are given for 4 stan- 
dard COMPTEL bands (0.75-1, 1-3, 3-10, 10-30 MeV). 
For a detection significance <2a the source flux is pre- 
sented by a 2a upper limit, otherwise by a flux point. 
Additionally, when comparing with the EGRET spec- 
tral extrapolations, the COMPTEL error bars and sub- 
sequently also the upper limits are enlarged by 20 percent 
to account for systematic errors. 

3. Results 

3.1. General Results 

We investigated all unidentified or tentatively identified 
EGRET sources for which spectral indices are given in the 
third EGRET catalogue. By excluding 5 artifacts near the 
Vela pulsar and 1 artifact near Crab l|Thompson et al.l 
l200lh . an ensemble of 173 sources was selected for our 
analyses. In some sky regions the unidentified EGRET 
sources are 'crowded'. For cases where they are closer 
than the COMPTEL location accuracy, no precise and 
unambiguous MeV flux could be attributed to individual 
EGRET sources. In these cases only one MeV source was 
fitted and its flux value was treated as upper limit for the 
different EGRET sources. 

In most cases, we derive only upper flux limits or marg- 
inal (<4cr) hints for MeV emission. We find significant 
MeV detections in 4 sky regions, three of them are loca- 
tions of already known unidentified MeV sources ((l,b)= 
(358.5, 0.5), (311.5, -2.5), (18.5, -0.5)). The fourth one, 
at ~(l,b)=(188.7, -4.4) and dubbed GRO J0550+19 in 
Bronsveld et al. (2002), is near the location of the Crab, 
which is by far the strongest MeV source. Detailed stud- 
ies of these sky regions have been reported ear l ier or 
are in progress llStrong et all EoOll IZhang et all l2002t 
iBronsveld et ~aTI. 120021: IC ollma.r et all l2004j) . 

In general our results allow the sources to be di- 
vided into 4 groups. For the majority of sources (~120) 
COMPTEL cannot provide any constraints on the 
EGRET spectra (group 1). We detect evidence (detection 
significance >2a) for about 20 sources whose MeV fluxes 
are consistent with the EGRET extrapolations, showing 
that the measured EGRET spectra extend into the MeV 
band (group 2). A third group contains a few sources 
indicating a spectral upturn of the EGRET spectra at 
MeV energies, suggesting the presence of an additional 
spectral component. Finally, we found 22 sources whose 
COMPTEL fluxes or flux upper limits are below the ex- 
pected fluxes based on extrapolations of their EGRET 
spectra, requiring a spectral turnover/break at MeV en- 
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Fig. 1. Combined simultaneous COMPTEL/EGRET en- 
ergy spectra of 4 sources out of our sample of 22 sources 
(see Table 1). The inserted text gives the EGRET 3rd cat- 
alogue informations (Hartman et al. 1999) on the source 
(source name, sky location, EGRET detection signifi- 
cance, and detection period). The solid line represents 
the best-fit EGRET spectrum above 100 MeV, the dashed 
lines its la error in spectral index, and the dotted lines the 
spectral extrapolations below 100 MeV down to 0.75 MeV. 
The required spectral changes at MeV energies are obvi- 
ous. The COMPTEL upper limits are 2a and the error 
bars on the flux points are la. 



ergies (group 4). In this paper we want to concentrate on 
the sources of group 4 for which COMPTEL can provide 
meaningful constraints on the EGRET spectra. Details on 
the three other groups will be given in a later paper. 

3.2. Sources with spectral constraints in the MeV band 

Applying the method described above, we found a sub- 
group of 22 unidentified (no tentatively identified source 
belongs to our sample) EGRET sources for which a spec- 
tral flattening at MeV energies is required. Fig. 1 shows 
four typical examples. All 22 sources are listed in Table 1 
including the relevant source parameters. At least two 
sources, 3EG J1638-5515 and 3EG J1823-1314, have a 7- 



ray luminosity peak in the 1-100 MeV band. Their lumi- 
nosity or luminosity upper limit in the COMPTEL 1-3 
MeV band is significantly lower than measured luminos- 
ity values at energies between 10 and 100 MeV, i.e. the 
COMPTEL spectral shape is harder than E" 2 , while the 
EGRET one is softer than E~ 2 . In fact, most likely the 
majority of sources in this sample has its maximum lu- 
minosity in this range, since we know that generally no 
strong hard X-ray sources have been found in EGRET 
error boxes. 

The spatial distribution of the 22 sources is shown in 
Fig. 2. A concentration at low galactic latitudes of |b| < 
30° is apparent. The sources tend to concentrate in the 
inner galactic region. Such a distribution suggests galactic 
origins for most of these sources. 

Fig. 3 shows correlation plots of spectral index ver- 
sus flux at energies above 100 MeV for a) the 22 selected 
sources and b) the rest (~151 sources) of the unidenti- 
fied and tentatively identified EGRET sources. No obvious 
correlation between spectral index and flux is visible for 
both source groups. The linear fit results in average pho- 
ton indices of 2.72+gjj for the 'break' sample and 2.13 ± 
0.03 for the rest. The softer energy spectra for the 'break' 
sample might be a selection effect due to the worse detec- 
tion sensitivity of COMPTEL compared to EGRET. For 
sources having similar flux levels in the EGRET band, 
those with a softer energy spectrum are more likely to 
be constrained by the COMPTEL data, because the ex- 
trapolation of the EGRET spectrum into the COMPTEL 
band will reach higher flux values. For the harder sources, 
the spectral extrapolation might go below the COMPTEL 
sensitivity limits making COMPTEL constraints impossi- 
ble. 



3.3. Variability 

There are three approaches to estimate the flux vari- 
ability of 7-ray source s. The so-called V method 
l|McLaughlin et all Il996(l is based on the x 2 -test. It is 
affected by the source detection significance, which is in 7- 
ray analyses/data often low. For a low-significance source 
with intrinsic variability a small V value is derived. On the 
other hand, a large V value can either be due to real source 
variability or due to large sys tematic effects. To over- 
come these systema tic problems, IZhang et all l|2000al) and 
iTorres et all Ij200ll) suggested to use the so-called I-index 
instead. It is defined as the ratio of the measured source 
variability to those of 7-ray pulsars, I=/i SO urcc/Mpuisars, 
which are considered to be intrinsically constant, fi is the 
ratio of the standard deviation of the measured flux values 
to the weighted mean flux. We regard a source as vari- 
able, if I-index > 2.5 (corresponding to a 3a significance 
for variability), a source as stable, if I-index < 1.5, and a 
source as dubious, if its I-index value is in between. We 
calculated t he I-indices for the 22 sources of our sample 
according to ITorres et alJ l)200lh by using the flux values 
given in the 3rd EGRET source catalogue. The results are 
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Fig. 2. The sky distribution of 173 unidentified and tentatively identified EGRET sources (+, A) in galactic coor- 
dinates. The triangles (A) represent the 22 EGRET sources of our sample. A clustering towards the inner galaxy is 
indicated. Especially, 20 of the 22 sources are located at |b| < 30°, suggesting a galactic origin. 



listed in Table 1. We found 10 sources to have an I-index 
> 2.5 (i.e. variable) and two to have an I-index < 1.5 (i.e. 
constant). 

iTompkinsI l)l999() investigated the source variability via 
the r-method, which is defined as r=cr//i, where a is the 
standard deviation of the fluxes and /i the average value. 
The idea is to overcome systematic effects by compar- 
ing to the t distribution for sources of known nature. 
The source fluxes are assumed to have a Gaussian dis- 
tribution with the parameters rx/i (width) and fi (mean 
value). These parameters are estimated by the maximum 
likelihood method. By comparing to the r distribution of 
blazars, a source is regarded variable if its 68% lower limit 
is greater than 1 . This results in four of our sources being 
variable. Three of them have an I value > 2.5 and one, 
3EG J1825+2854, has an I value of 2.49. These results 
show that the I- and r-methods provide consistent results. 
All other 18 sources, with r max values larger than 0.3 and 
T m in values less than 1, wo uld be class i fied a s dubious, ac- 
cording to the definition bv lTompkinsI (^999). Nolan et al. 
(2003) adopted the same r-method, only applying slight 
changes/modifications, to investigate the time variability 
of the EGRET 7-ray sources. Their approach resulted 
in some numerical differences to Tompkins. According 
to the classification by Nolan et al. (2003), four sources 
of our Table 1, 3EG J1612-2618, 1828+0142, 1837-0423, 
1940-0121 and 2034-3110, are variable. These sources have 
I-indices > 2.5. All others, with the exception of 3EG 
J1424+3734, which is not included in Nolan et al. (2003) 
due to their selection criteria on the observations, are 
classified as dubious. Since the t values are estimated 
by lumping the VPs within one month while the I-index 
method uses individual VPs, typically ranging from sev- 
eral days to two weeks, the two methods may estimate 
the source variability on different time scales. We use the 
I-index to classify the source variability. This yields, as 



mentioned above, 10 variable sources. Because of the large 
I-index value (2.49) and the variability-indicating r value 
(^min ~2.59), we consider also 3EG J1825+2854 as a vari- 
able source. The two sources with I-index < 1.5 are consid- 
ered dubious due to their r values. Thus we have a sample 
of 11 variable sources and of 11 dubious sources. 

Fig. 4 shows correlation plots of spectral photon index 
versus I-index for the 22 sources of our sample and the 
remaining 151 sources. For the sample of the 151 sources, 
the linear correlation coefficient is derived to be ~0.15, 
suggesting that the softer sources tend to be more variable. 
This is consistent with Torres et al. (2001), who found that 
the most-variable unidentified EGRET sources near the 
galactic plane tend to have steep spectra. The 22 sources 
do not show such a trend. A linear correlation coefficient of 
~0.03 indicates an uncorrelated sample. This might either 
discriminate this sample of 22 sources from the rest, or is 
due to the relatively poor statistics. 

4. Discussion 

By analyzing the COMPTEL data of unidentified and ten- 
tatively identified EGRET sources, and comparing their 
COMPTEL spectra with those from EGRET, we found a 
subgroup of 22 sources which show a spectral turnover or 
break at MeV energies. At least two of them, but most 
likely the majority, have an emission maximum at ener- 
gies between 1-100 MeV. Variability analyses show that 
half of the sources are variable above 100 MeV. For the 
other half no significant time variability could be proven. 

The properties of the variable sources are reminis- 
cent to the observational properties of EGRET-detected 
blazars. They are generally variable and their broadband 
(radio to 7-ray) spectra are characterized by two spec- 
tral maxima. The high-energy one is for many sources lo - 
cated at MeV energies, e.g. 3C 273 l(Lichti et all Il995h . 
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Table 1. Catalogue of the 22 unidentified EGRET sources which - by inclusion of the COMPTEL results - have 
to have a spectral break or their 7-ray luminosity maxima at MeV energies. The source parameters of the first six 
columns are taken from the third EGRET catalogue, column 7 and 8 are for the variability analysis, and the last 4 
columns give the COMPTEL fluxes/upper limits in the 4 standard energy bands. The abbreviations are the following: 
Name - 3EG source name. 

I. b. - source coordinate of galactic longitude and latitude. 

F - EGRET source flux and error in units of 10 -8 ph cm~ 2 s _1 . 

a - photon spectral index and its error. 

VP - CGRO observation periods for which the spectral index is derived. PI means phase 1, P12 the combination of 
PI and P2, and so on. Ob servation 330+ is the combined VPs 330 and 332. 

t (Tmin ,Tmax ) - taken from iTonrnkinsI l)l999() for measuring the fractional variability of EGRET sources. T mm and r max 
represent separately the 68% low e r and upper limits on r. 

/ - index defined bv lTorres et"aD l)200l|) for measuring source variability with respect to 7-ray pulsars; derived by our 
variability analysis on the published EGRET fluxes. 

Flux (10~ 5 ph cmT 2 s' 1 ) - COMPTEL fluxes and the flux unit in the four standard energy bands: 0.75-1, 1-3, 3-10 
and 10-30 MeV. The upper limits are 2a and the error bars are la. 



Name 


1 


b 


F 




a 




VP 


T ( 7~m i n 5 T~m ax ) 


I 


Flux (10 -5 ph 


cm J s 


L ) 


3EG J 


deg 


deg 


















0.75-1 


1-3 


3-10 


10-30 


0429+0337 


191.44 


-29.08 


12.0+2.7 


3 


02+0 


27 


P1234 


0.00(0.00,0.45) 


1 


6 


<3.2 


<4.0 


<3.1 


<1.2 


0520+2556 


179.65 


-6.40 


15.7+2.7 


2 


83+0 


24 


P1234 


0.00(0.00,0.31) 


1 


1 


4.8+1.7 


8.3+1.6 


<2.1 


<0.5 


0546+3948 


170.75 


5.74 


13.7+2.6 


2 


85+0 


21 


P1234 


0.11(0.00,0.47) 


1 


7 


<2.9 


<2.9 


<1.4 


<0.9 


1300-4406 


304.6 


18.74 


10.6+2.9 


3 


07+0 


40 


P12 


0.48(0.00,1.57) 


3 





<5.0 


<4.7 


<2.2 


<1.0 


1424+3734 


66.82 


67.76 


16.3+4.9 


3 


25+0 


46 


PI 


0.01(0.00,oo) 


1 


9 


<6.4 


8.5+3.1 


<5.6 


<2.2 


1500-3509 


330.91 


20.45 


10.9+2.8 


2 


99+0 


37 


P1234 


0.00(0.00,0.61) 


1 


5 


<5.1 


<3.3 


<1.6 


<1.0 


1612-2618 


349.40 


17.90 


92.2+27.7 


2 


71+0 


23 


423. 


1.78(0.76,11.74) 


4 


1 


<41.0 


<19.1 


<9.3 


<4.2 


1638-5515 


334.05 


-3.34 


67.3+14.2 


2 


56+0 


21 


P2 


0.00(0.00,0.69) 


2 


4 


<15.6 


<7.0 


<6.7 


2.8+1.0 


1639-4702 


337.75 


-0.15 


53.2+8.7 


2 


50+0 


18 


P1234 


0.00(0.00,0.38) 


2 





<3.1 


<4.5 


<3.0 


<0.9 


1709-0828 


12.86 


18.25 


12.6+3.2 


3 


00+0 


35 


P1234 


0.84(0.11,2.21) 


2 


7 


<3.2 


<3.2 


<2.8 


<0.6 


1735-1500 


10.73 


9.22 


196.3+48.8 


3 


24+0 


47 


231.0 


1.09(0.00,10.14) 


8 


9 


<29.0 


<23.1 


<11.8 


<4.5 


1741-2312 


4.42 


3.76 


33.1+5.9 


2 


49+0 


14 


P12 


0.52(0.18,1.03) 


2 


2 


<11.9 


<4.4 


<3.0 


<0.9 


1800-0146 


25.49 


10.39 


26.1+6.1 


2 


79+0 


22 


P34 


0.00(0.00,0.48) 


1 


9 


<8.1 


<7.0 


<2.2 


<1.4 


1823-1314 


17.94 


0.14 


102.6+12.5 


2 


69+0 


19 


P3 


0.72(0.40,1.37) 


3 





10.2+2.9 


5.6+2.8 


3.7+1.2 


2.7+0.5 


1825+2854 


56.79 


18.03 


34.3+10.9 


4 


47+1 


15 


9.2 


73.19(2.59,oo) 


2 


5 


<30.1 


<13.9 


<8.6 


<2.2 


1828+0142 


31.90 


5.78 


132.2+24.0 


2 


76+0 


39 


13.1 


3982.15(6.92,oo) 


5 


3 


<16.4 


<15.8 


<11.3 


<4.4 


1837-0423 


27.44 


1.06 


310.4+63.7 


2 


71+0 


44 


423.0 


12.01(2.17,oo) 


8 


4 


<23.4 


<18.5 


<14.3 


<3.0 


1858-2137 


14.21 


-11.15 


11.2+2.6 


3 


45+0 


38 


P1234 


0.00(0.00,0.56) 


2 


8 


<4.2 


<5.0 


1.6+0.7 


<0.6 


1903+0550 


39.52 


-0.05 


62.1+8.9 


2 


38+0 


17 


P1234 


0.35(0.18,0.60) 


2 


3 


<3.3 


<3.4 


<2.7 


<0.7 


1940-0121 


37.41 


-11.62 


41.0+10.7 


3 


15+0 


39 


330.+ 


4.58(1. 13,oo) 


4 





<16.6 


<10.9 


<6.8 


<1.8 


2020-1545 


28.09 


-26.62 


11.8+3.4 


3 


40+0 


55 


PI 


0.00(0.00,0.80) 





9 


<5.2 


<7.7 


<2.6 


<1.0 


2034-3110 


12.25 


-34.64 


17.4+5.2 


3 


43+0 


78 


PI 


2.88(0.89,154.84) 


5 


7 


<5.4 


<5.4 


<3.2 


<1.7 



PKS 0528+13 4 llCollmar et al l ll 997t iMukneriee et~afl. 
Il999l) . 3C 279 l|Hartman et all 1200 ill. Fo r detai ls on their 
high-energy emission see e.g. ICollmarl ()200l[l and ref- 
erences therein. Because the EGRET blazars are typ- 
ically high-latitude sources, variable unidentified high- 
latitude EGRET sources are generally considered to be 
of extra-galactic blazar origin. Due to EGRET's de- 
creasing sensitivity to variability with decreasing flux 
they can appear as non-v ariable sources at low latitudes 
( McLaughlin et all flQQfih . where, due to the strong dif- 
fuse galactic emission, a higher significance level for source 
detection is required. Therefore some of the low-latitude 
unidentified EGRET sources may also be blazars. Two 
sources of our sample, 3EG J1424+3734 and 3EG J2034- 
3110, are high-latitude (|b| > 30°) sources. 3EG J2034- 



3110 is variable and 3EG J1424+3734 is dubiously vari- 
able. Our finding of a 'blazar-like' MeV spectrum provides 
further evidence for the blazar nature of these 2 sources, 
especially for the case of 3EG J2034-3110, for which time- 
variability is observed. The lack of an obvious blazar can- 
didate at thes e source posi t ions c an be explained by the 
calculations of lTorres et al . ( 2002). The unidentified 7-ray 
sources could be a distant weak 7-ray emitting blazars 
whose emission is amplified by gravitational microlcnsing. 
The magnification factor is energy dependent and a bro- 
ken energy spectrum is predicted. Although the individual 
source properties of our sample sources are reminiscent of 
blazars, their distribution on the sky is not. While the 
blazars are mainly detected at high latitudes, our source 
group concentrates towards the inner galaxy. Therefore we 
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Fig. 3. Correlation plots of spectral photon index versus 
flux (EGRET results) for our sample of 22 sources (upper 
panel) and the rest of 151 unidentified and tentatively 
identified EGRET sources (lower panel). The sources of 
our sample typically have softer spectra. 



conclude, that at least the majority of these sources are 
not of blazar origin. 

Some 7-ray pulsars are also showing such an MeV- 
turnover spectrum . One example is PSR B1509-58. 
iKuip er et alJ l)l999j) showed that its maximal luminosity 
is reached between 10 and 30 MeV, followed by a strong- 
softening of its spectrum towards the EGRET band. In 
fact, PSR B1509-58 has not been detected by EGRET 
above 100 MeV, while it is a strong COMPTEL source. 
PSR 1509-58 is a young pulsar and its surface magnetic 
field can be inferred as at least 3.1xl0 13 G (Kuiper et 
al. 1999). In such sources, electrons are accelerated up to 
relativistic energies, and subsequently produce 7-rays via 
curvature radiation and inverse-Comptonization of soft 
photons. In the vicinity of the magnetic pole, 7-rays can 
be absorbed by the strong magnetic field via the pho- 
ton splitting process. The latter process happens if the 
magnetic field is larger than a critical value B > 0.3 B cr , 
where B rr =4.413x 10 13 G is the surface magnetic field 
ijHarding et all fl997j) . Its onset has no energy threshold. 
For such a strong magnetic field, the polar cap scenario 
predicts that photon splitting can become the dominant 
attenuation process, resulting in a softening of the energy 
spectrum below 100 MeV, i.e. providing the observed MeV 
cutoff. It should be mentioned that also in the competing 
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Fig. 4. Correlation plots of spectral photon index ver- 
sus variability I-index for our sample of 22 sources (upper 
panel) and the rest of 151 unidentified and tentatively 
identified EGRET sources (lower panel) . The lower frame 
indicates a trend for higher variability with softer spectra. 

outer gap scenario such a spectral break at low 7-ray en- 
ergies can be accounted for. Namely, Zhang and Cheng 
(2000b) applied their three dimensional outer magneto- 
sphere model to PSR B1509-58. They could reproduce the 
measured broad pulse profile and the measured pulsed en- 
ergy spectrum from the optical range up to 7-rays. This 
spectrum shows a power-law shape with a spectral bend 
above 1 MeV. In the same work, Zhang and Cheng consid- 
ered also the case of the young Crab-like LMC pulsar PSR 
B0540-69, which is also not detected by EGRET but below 
10 keV stronger than PSR B1509-58. The overall charac- 
teristics of these two young pulsars appear rather similar. 
By analyzing RXTE data, De Plaa et al. (2003) showed 
that PSR B0540-69 and the Crab pulsar have very similar 
spectral shapes up to about 50 keV, where the Crab pul- 
sar spectrum reaches its maximum luminosity. From these 
hard X-rays to the EGRET energies above 100 MeV, the 
Crab spectrum softens throughout the COMPTEL range 
(photon index ~ -2.4; Kuiper et al. 2001). De Plaa et al. 
discussed the spectral shapes of these three very young 
pulsars (< 1.6 x 10 3 yr) which appear to be different 
from those of older 7-ray pulsars like Vela and Geminga. 
The youngest pulsars are strong(er) in the X-ray domain, 
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but weak(er) above 100 MeV. Such a spectral behaviour 
makes them less likely counterparts to the unidentified 
EGRET sources, because their Lx/L 7 ratio is too high. 
On the other hand, all older established 7-ray pulsars 
like Vela, PSR B1706-44, PSR B1951+32, Geminga and 
PSR B1055-52 are weak X-ray emitters but have harder 
7-ray spectra extending up to the GeV range. The lat- 
ter might be a selection effect. For example, the weakly 
detected (above 50 MeV) 7-ray pulsar PSR B0656+14 ex- 
hibits a steep s pectrum (index = -2 8±0.3 ) at energies 
above 100 MeV (Ramanamur This pulsar 

has a modest surface magnetic field of 4.7 x 10 12 G and a 
characteristic age of 1.1 x 10 5 years. It has not been de- 
tected by COMPTEL below 30 MeV, implying a turnover 
of the soft EGRET spectrum before the COMPTEL win- 
dow, just like our sample. 

Another pulsar type, the old, recycled weakly mag- 
netized millisecond pulsars, can also show such an MeV- 
turnover spectrum. Namely, a 4.9a detection has been 
claimed of the 2.3 ms pulsar PSR J0218+4232 (Kuiper et 
al., 2000; 2002) at EGRET energies. Its spin down param- 
eters give a rather old characteristic age of 4.6 x 10 s years. 
Its magnetic field of 4.3 x 10 s G is much weaker than that 
of standard 7-ray pulsars (e.g. Crab, Geminga). A weak, 
but very hard X-ray spectrum was measured (photon in- 
dex below 10 keV > -1; Minco et al. 2000) and a rather 
soft 7-ray spectrum (photon index ~ -2.6) was detected for 
EGRET above 100 MeV. Current polar cap and outer gap 
models can account for the production of such high-energy 
radiation in millisecond pulsars, but have difficulties in re- 
producing the observed spectral shape. Nevertheless, the 
observed spectrum has the overall spectral shape of our 
source sample, making millisecond pulsars viable candi- 
date counterparts, but only for the stable EGRET sources. 
Millisecond pulsars are expected to be stable high-energy 
emitters just like the normal radio pulsars, and contrary to 
case of accreting pulsars for which variability is naturally 
expected. 

X-ray binaries (XRBs) have also been suggested as 
counterparts of unidentified EGRET sources. One XRB, 
the neutron star system Cen X-3, has been detected as 
7-ray emitter during an activity period in 1994 (Vestrand 
et al. 1997). EGRET found a temporary 5<r source which 
was positionally coincident with Cen X-3. The 4.8s mod- 
ulation of the 7-rays, coinciding with the 4.8s rotation pe- 
riod of the neutron star, provided compelling evidence for 
the identification (Vestrand et al. 1997). A hard power- 
law spectrum (photon index -1.81±0.37) was measured 
between 70 MeV and 10 GeV. The authors suggest that 
galactic X-ray binary systems may constitute a class of 
highly variable GeV 7-ray sources. It is assumed that in 
Cen X-3 we see the unabsorbed 7-ray spectrum. If in such 
a system the radiation region is surrounded by a con- 
densed soft photon field, e.g. in an accreting XRB, an 
energy-dependent absorption of the 7-rays will occur re- 
sulting i n a soft 7-ray spectr um. This scenario was pro- 
posed bv lRomero et al.l l 200ll) for the possible association 
of the EGRET source 3EG J0542+2610 and the Be/X-ray 



transient A0535+26. Such a n absorption process was also 
studied bv lWu et all ill993h for the XRB Cyg X-3. Their 
simulations showed that the 7-ray spectrum can change 
significantly when passing through the ambient soft X-ray 
field of an accreting source. The 100 MeV - 1 GeV emis- 
sion will be absorbed. For both cases a soft 7-ray spectrum 
would be observable by EGRET, if the absorption is not 
too strong. 

The emission processes of microquasars/blazars and 
extragalactic blazars are - in principle - the same, how- 
ever, on different time, space, and energy scales. Therefore 
one expects microquasars as pote ntial counterp a rts o f 
the unidentified EGRET sources. IParedes et ail l)2000h 
suggest the microquasar LS 5039 to be the counterpart 
of 3EG J1824-1514. Calculations of microquasar spectra 
show that, depending on the strength of the jet Lorentz 
factor and magnetic field, microquasars could be detected 
by EGRET, and that they could have their spectrum 
turnover at MeV energies (e.g. Kaufman Bernadd et al. 
2002). However, no microquasar was definitely identified 
yet as an EGRET source. Given the current knowledge, 
the properties of some of our sources, spectral turnover 
and variability, could be matched by microquasars. 

Binaries composed of early-type stars, like Wolf Rayet, 
Of and Be stars, which produce strong stellar winds are 
proposed by Benaglia & Romero (2003) to be potential 
7-ray emitters. They argued that the electrons could be 
accelerated to relativistic energies by the shocks gener- 
ated in the colliding wind region of the early-type bina- 
ries, and then cool via the process of inverse-Compton 
scattering off the local soft photon field. A low-energy cut 
off of the electron spectrum or incomplete cooling of the 
electron population could lead to a turn over in the 7-ray 
spectrum at lower 7-ray energies. The early-type binaries 
are concentrated in the inner spiral arms of the Galaxy 
and variability is expected due to the changing geometry. 
However, according to the investigation of Romero et al. 
(1999) on the spatial correlation of low-latitude uniden- 
tified EGRET sources with early-type stars, none of the 
22 sources of our sample is spatially coincident with any 
early-type binary. 

Another class of sources proposed to be candidates 
for unidentified EGRET sources are supernova remnants 
(SNRs) (e.g. Esposito et al. 1996). In SNRs electrons 
and/or protons can be accelerated by the Fermi process 
to relativistic energies or appear in the outflows of a pul- 
sar, if one is embedded in the SNR. Subsequently 7-rays 
can be produced via processes of inverse-Compton scatter- 
ing, relativistic bremsstrahlung, synchrotron emission, or 
7To decay. Since the first-order Fermi process has a mono- 
tonic evolution of energies with time and the acceleration 
time for the maximum energy is limited by the age of 
the remnant, it is generally thought that SNRs should 
be stable 7-ray sources, at least on a time scale of sev- 
eral years. SNRs arc spatially large, have ambient sur- 
rounding matter and a relatively weak magnetic field. A 
spectral turnover at MeV energies as found for our sam- 
ple could be generated by inverse-Comptonization of soft 
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photons, by having the maximum of the synchrotron emis- 
sion at MeV energies (e.g. Crab nebula), or by the decay of 
7To particles (generated in p-p collisions). The latter one 
would result in a broad spectral bump centered around 
68 MeV. Four sources of our sample, 3EG J1639-4702, 
1823-1314, 1837-0423, 1903+0550, are positionally coinci- 
dent with SNRs. 3EG J1823-1314 and 3EG J1837-0423 are 
variable (Table 1) and therefore are unlikely counterparts 
of the SNRs. For 3EG J1639-4702 and 3EG J1903+0550 
no variability is proven, therefore they remain potential 
counterparts for the EGRET 7-ray sources. Their possi- 
ble associations with SNRs and coinciding radio pulsars 
have been discussed by Torres et al. (2003). With respect 
to radio pulsars, they concluded based on the energetics, 
that for 3EG J 1639-4702 only one of the coinciding radio 
pulsar, PSR J1637-4642, has the potential to be the coun- 
terpart. For 3EG J1903+055 they did not find a potential 
pulsar counterpart, although they note that one of the co- 
inciding pulsar lacks information and therefore can not be 
judged. With respect to 7-ray emission due to 7ro-decay, 
Torres et al. point out that the SNRs are too far away to 
generate the observed flux. Even considering an enhance- 
ment of the 7-ray production due to interactions of the ac- 
celerated nuclei with nearby dense molecular clouds, they 
find the identifications to be unlikely. These two unident- 
fied sources could be associated with a pulsar wind nebula, 
analogous to the case of the Crab for which the electrons 
might be accelerated in the inner nebula and subsequently 
are generating a synchrotron spectrum which cuts off at 
MeV energies (de Jager et al.,1996). Their fluxes above 
100 MeV have about the level of the Crab nebula. Because 
with distances of at least 8 kpc (Torres et al., 2003), i.e. at 
least 4 times the distance to the Crab, they would have to 
be significantly more powerful at EGRET energies than 
the Crab nebula. However, this could be possible, if their 
synchrotron cut off would be shifted to higher energies 
(higher characteristic synchrotron energy) compared to 
the Crab nebula. 

Because for many unidentified EG RET sources no ob- 
vious candidate counterpart exists, iRomero et alJ l)l999(l 
mentioned isolated rotating black holes, standard (Kerr) 
or charged (Kerr-Newman) ones, as possible counter- 
parts. They could accrete from the diffuse interstel- 
lar medium (ISM). Changes in the density of the ISM 
would result in a variable 7-ray flu x. One such s ource , 
3EG J1828+0142, was modeled by IPunslv et all l|2000|) 
by assuming an isolated Kerr-Newman black hole ori- 
gin. The model predicts a steep synchrotron self-Compton 
spectrum with a spectral maximum at MeV energies, i.e. 
matching the spectral properties of our source sample. 
Including 3EG J1828+0142, seven of our sample sources 
(3EG J0520+2556, 0546+3948, 1638-5515, 1735-1500, 
1741-2312, 1800-0146, 1828+0142) are located at low 
galactic latitudes (|b|< 10°) and lack any positional coin- 
cidence with galactic objects of potential counterpart na- 
ture. Two of them, 3EG J1735-1500 and 3EG 1828+0142, 
are significantly variable and therefore match the antici- 
pated 7-ray properties of these exotic objects. 



5. Conclusion 

By analyzing the contemporary COMPTEL observations 
of all unidentified or only tentatively identified EGRET 
sources, we found a subgroup of 22 sources for which 
we can provide spectral constraints for source modelling. 
Their spectra have to turn over between ~1 MeV and 
100 MeV, and at least two of them, but most likely the 
majority, have their maximum luminosities somewhere in 
this energy band. Most of the sources are not detected 
by COMPTEL, however the simultaneously derived up- 
per limits require the spectral bending. These sources have 
rather steep energy spectra in the EGRET band, and seem 
to be preferentially located in the inner galaxy, especially 
at low latitudes (|b| < 30°). Variability studies reveal that 
half of them are significantly variable above 100 MeV. 
Potential counterparts have to conform to these observa- 
tional results. A blazar origin for the two high-latitude 
sources in this sample seems to be likely. Viable candidate 
counterparts for the steady low-latitude sources are: 1) 
young (age < 10 6 years) pulsars, although the youngest 
with the strongest magnetic fields might be too strong 
in the X-ray domain, 2) old, recycled millisecond pulsars 
with a weak magnetic field (like PSR J0218+4232), and 3) 
SNRs and pulsar wind nebula whose synchrotron spectra 
are peaking at MeV energies (like the Crab nebula). For 
the variable low-latitude sources, XRBs, in particular mi- 
croquasars/blazars by assuming a spectral analogy to the 
extragalactic objects, and isolated BHs would match the 
requirements. Case by case studies might reveal further 
insights in the nature of individual sources. 
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